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Ergodicity of the systems with Nose-Hoover thermostat are studied. The dynamics of the heatbath variables are investigated and they 
can be periodic when the system has quick oscillation. The periodic behaviour of them causes the system to lose its ergodicity. The 
kinetic-moments method is also studied, and the heatbath variables in this method are found to be chaotic. The chaotic behaviour makes 
the whole system ergodic. 
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1 Introduction 



| Molecular dynamics (MD) simulations have become more powerful tools with the development in per- 
■ formance of computers. The standard MD simulations have been performed with constant energy, and 
*~Q I therefore, these simulations have taken samples in the microcanonical ensemble. It is required to per- 
form constant-temperature MD (Q), because the most of physical quantities are observed in constant- 
ly ■ temperature environment. Nose proposed the extended system method which allows to simulate constant- 
| temperature MD conserving time-reversible property (J3). Hoover reformulated the Nose's method to a 
more practical one, which is now called the Nose- Hoover method 
Recently, the Nose-Hoover thermostat has been applied not only for the equilibrium system but also for 
' non-equilibrium phenomena. Ogushi et al. studied the heat conduction of the three-dimensional Lennard- 
Jones particle system using molecular dynamics simulation (jj). They put two thermostat with different 
temperatures on the both ends of the system to achieve the spontaneous phase separation. 
The Nose-Hoover method achieves Gibbs' canonical ensemble when the system is ergodic. However, some 
O ■ numerical researches report that the system with a small degree of freedom often loses its ergodicity, and 
consequently, the system cannot be canonical in Gibbs' sense. In the system studied by Ogushi et al, for 
instance, the density at the higher-temperature side becomes dilute, and the number of particles controlled 
^ ', by the thermostat can be very small. Therefore, the system can lose the ergodicity. 

In order to improve the ergodicity of a system with thermostat, many methods have been proposed 
However, there are less studies investigating why and how the system loses its ergodicity. In the present 
article, we focus the ergodicity of some systems in which temperature controlled by the Nose-Hoover 
method. 
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2 Nose-Hoover Method 

Consider an isolated system described by a Hamiltonian 7~C(p, q) with momenta p and coordinates q. The 
equations of motion are 

■ m m 

Qi = "5— > (!) 

= — «-i (2) 
% 

where ^4 denotes the differentiation of A with respect to time. The energy of the system is conserved. 

Following Nose and Hoover, we can achieve the canonical ensemble by modifying the equations of motion 
as follows: 



dU 

dpi 



(3) 



Pi = -o CPi, 4 ) 

r, = (, (6) 

where £ and r/ are additional variables expressing the Nose-Hoover thermostat, and K, Kq, and r are the 
kinetic energy, the aimed temperature, and the relaxation time of the thermostat, respectively. 
While the system with the Nose-Hoover thermostat does not have the canonical form, a quantity 

H' = H + K (r( 2 + 2 V ) (7) 

is conserved with the original Hamiltonian TL. The energy of the system without the thermostat can 
fluctuate with the Boltzmann weight exp(— (3TL). Therefore, the system will achieve the canonical ensemble 
when the system is ergodic. 



3 Simulation 

In order to investigate the ergodicity of the system with the Nose-Hoover thermostat, we study the following 
two systems, 

^ho = ^ 2 + 9 2 ), ( 8 ) 
n y ^ = \(p\+pl) + \qlql (9) 

where T^o an< ^ ^y m are ^ ne Hamiltonian of the harmonic oscillator and the Yang-Mills type, respectively. 

We performed numerical simulations of the two systems with the Nose- Hoover thermostat. The time 
evolution of the systems were calculated by the fourth-order Runge-Kutta method. The time step At is 
set to be 0.005, the aimed temperature Kq is 0.5, and the total step is 10 7 . The results for the harmonic 
oscillator are shown in Fig. This figure shows that the system loses the ergodicity, and the distribution 
of the energy is far from that of the canonical ensemble. 

The results of the Yang- Mills system are shown in Fig. |2j The phase diagram of (pi,P2) is chaotic, and 
the histogram of the energy shows canonical distribution, while the harmonic-oscillator system does not. 
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Figure 1. (a) Phase diagram and (b) an energy histogram of the harmonic-oscillator system. 




Figure 2. (a) Phase diagram and (b) an energy histogram of the Yang-Mills system. 

It suggests that the Yang-Mills system with the Nose-Hoover thermostat is ergodic while it has only two 
degrees of freedom. 

We focus the dynamics of the heatbath variables in order to figure out why the two systems show different 
behaviours. The time evolution of the heatbath variables for the harmonic-oscillator system is shown in 
Fig. 01 The both variables seems to be periodic, and the absolute values of them are bounded throughout 
the simulations. 

The time evolution for the Yang-Mills system is shown in Fig. 0] The behaviour of the heatbath variables 
seems to be chaotic, and the amplitude is much larger compared to the harmonic-oscillator case. The 
chaotic behaviour of the heatbath variables is essential for the ergodicity of the system with Nose-Hoover 
thermostat. 



4 Dynamics of Heatbath variables 

In this section, we derive the differential equation for the heatbath variables in order to investigate why the 
variable show the periodic behaviour in the harmonic-oscillator system. The equations of motion Eqs.(|5j). 
(@J), and © are transformed to be 

c+2cc+4k = -^^, (io) 

in the harmonic-oscillator case. 
With large enough t, the heatbath variable £ varies much slower than pq. Therefore, we can take an 
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Figure 3. The harmonic-oscillator system. Time evolutions of the heatbath variables (a) £ and (b) rj are shown. 




Figure 4. The Yang-Mills system. Time evolutions of the heatbath variables (a) £ and (b) r] are shown. 

approximation {pq) = 0. Assuming the scaling form for £ to be 

at) = -at/r), (ii) 

r 

we have the differential equation for £ to be 

C + 2CC + C = 0, (12) 

which is independent of r. Therefore, r-dependence of £ should have the form in Eq. (jllj) . A quantity 

#C = C 2 + <Mog(^+l) (13) 

is conserved in Eq. (|12|h We can prove that £ and £ are bounded because of £ + 1 = K/Kq > 0. Therefore, 
the variable Q should be periodic. Note that r\ should also be periodic since C is an odd function. 

We solved Eq. (|12|) numerically using the fourth-order Runge-Kutta method. The initial condition is 
C\t=o = 0, C|t=o = 1- The results are shown in Fig. The behaviour is periodic as expected. 

5 Kinetic-Moments Method 

Kinetic-moments method is one of the methods proposed to improve the ergodicity even for the harmonic- 
oscillator case (ja). This method controls the kinetic energy (K) and its fluctuation (-f^ 2 ) — {K) 2 simulta- 
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Figure 5. Time evolution of f. The initial condition is C|t=o = 0, C|t=0 = 1- The behaviour of the heatbath variable is periodic. 
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Figure 6. Time evolutions of the heatbath variables, (a) Time evolution of £ and (b) that of £ are shown. The behaviour seem to be 
chaotic and the amplitudes of the fluctuation are much larger than that of the Nose-Hoover system. 

neously. The equations of motion in this method are 



P 



on 

dp ' 
&H 
dq 



( = {p 2 - i) /r c , 

i = (p 4 - 3p 2 ) /r € , 



(14) 

(15) 

(16) 
(17) 



with an additional heatbath variable £. While Hoover and Holian reported the system with the kinetic- 
moments thermostat shows good canonical distribution, they did not explain why the method achieves 
ergodicity. 

Here, we discuss the behaviour of the heatbath variables in the method. We simulated the harmonic- 
oscillator system with the kinetic-moments thermostat. The details of the simulation are same as the 
previous. The time evolution of the two heatbath variables £ and £ are shown in Fig. H3 The heatbath 
variables seem to be chaotic while the variables are periodic in the Nose-Hoover system. 

We can derive the differential equations for the heatbath variables in the kinetic-moments case. The 
equations are complex, and are expected not to have any conserved values. Therefore, the kinetic-moments 
method achieves the ergodicity because of the chaotic behaviour of the heatbath variables, not by control- 
ling higher moments of the kinetic energy. 
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6 REFERENCES 
6 Summary and Discussion 

We have studied the ergodicity of the Nose-Hoover thermostat in the two systems, the one-dimensional 
harmonic-oscillator and Yang-Mills systems. The numerical results have shown that the Yang-Mills system 
is ergodic while the harmonic-oscillator system is not. 

In order to study why the two systems have shown the different behaviour, we have derived the differential 
equation for the heatbath variable. When the system oscillates quickly, the heatbath variable becomes 
periodic, and consequently the whole system cannot be ergodic. Therefore the Nose-Hoover system loses 
its ergodicity not because the system has a small degrees of freedom, but the heatbath variables becomes 
periodic. A particle system with the Nose-Hoover thermostat will be ergodic, since quick oscillations are 
not expected to exist, even if the system has small degrees of freedom. We also explain r-dependence of 
the heatbath variable in the amplitude and the period at the same time. 

We have also studied the ergodicity of the kinetic-moments method. This method contains an additional 
heatbath variables to the Nose-Hoover method, and it makes the heatbath variables chaotic. The chaotic 
behaviour of the variables achieves the ergodicity of the whole system. These results imply that at least 
one of the original system and the heatbath variables should be chaotic in order that the whole system is 
ergodic. 
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